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A Reinvestigation, of the Problem of the Automorphic 
Linear Transformation of a Bipartite Quadric. 

By Thomas Muir, LL. D. 
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the problem is to find two matrices 
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we shall have 
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2. The mere performance of the substitutions changes the given bipartite 
quadric into 
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which by a property of bipartite functions,* 
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Consequently if we denote the square array of the original quadric by D, and 



* Trans. Roy. Soc. Edinburgh, XXXII, pp. 461-483. 
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those of the two substitutions by S lt /S 2 respectively, the transformed function 
may be written* 

£ n K 



-Si 



trD 



tr So 



e 



But this by another property of bipartites is equal to 

£ 37 z ... . 



tr So.D.S, 



e 



so that the problem is changed into finding two matrices S t and S 2 which will 

satisfy the equation 

trS t .D.S 1 = D. 

The problem is evidently indeterminate, there being two unknown matrices and 
one matrical equation to be satisfied, or 2n 2 scalar unknowns and n 2 scalar equa- 
tions for obtaining them ; the solution therefore should involve n 2 arbitrary 
scalar constants. 

3. It is immediately clear that we shall obtain a solution — and the most 
natural one — by putting S l or S 2 equal to any arbitrary matrix and then solving 
for S. z or iSy. Doing this and denoting the arbitrary matrix by M we have the 
alternative solutions 

Sy = M, I 

trS l = DM- 1 D-\) 

trB, = M, I 

S 1 = D- l M- 1 D.l 



*" tr D" is used, as by Cayley, for the ' transverse' of D, the symbol tr indicating the operation of 
row-and-column exchange: " conj. D" would have been more appropriate, considering the way in 
which ' conjugate ' is already used in connection with the elements, terms and minors of a determinant. 
29 
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4. A symmetry more in keeping with the character of the equation will be 
attained by putting 8 X equal to the product of an arbitrary matrix and D, and 
then as before solving for S 9 . This course of procedure gives us 

S 1 = MD, 

trS, 



or 



\ = MD, \ 

S 1 = M- 1 D,\ 
trS i = MD~\\ 



In all these forms of solution, it will be observed, M, though in other 
respects arbitrary, must not have a vanishing determinant. 

5. Taking one of the forms, say the third, we have the following simple 
theorem for the transformation in question : 

The transformation of 

(x, y,z \D\x', y', z' ) into (£, n , f , V>W, *?', £' ) 

can be effected by the substitutions 

(as, y, a,... .)= {MD l£ t *i,£, ), 

(x', y', *,... .) = tr{M-'D-^, ti', ? ), 

where Mis any arbitrary matrix ivith a non-vanishing determinant. 

6. For the case of an ordinary quadric, viz. where D is axisymmetric, equal 
to A say, and aS' 3 = S 1} = S say, the equation for solution is 

trS.A.S=A. 

Now the left-hand member, as well as the right-hand member, of this is axisym- 
metric; for by the law of row-and-column exchange in a product 

tr (tr S.A.S) = tr S. trA.tr (tr S) 
= trS.A.S. 

The state of matters here therefore is that though there is one unknown matrix 
and one matrical equation to determine it, there are n z scalar unknowns, and 
only I n (n + 1) scalar equations to be satisfied. The problem is thus again an 
indeterminate one, there being bn(n — 1) arbitrary constants to be expected in 
the solution. 

One solution of the equation is evident, viz. 

trS= — l: 
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it follows therefore that tr S + 1 is a factor of trS.A.S — A. Denoting the 
other factors — the first terms of which must evidently be A and S — by 

A + x and S + y , 

and bearing in mind that the product of the three is axisymmetric, we see from 
the law of row-and-column exchange in a product that 

y = l 

and x is axisymmetric. 

The first and third factors being thus known, the second is readily found to be 

A — A (1 - S)- 1 — (1 — tr S)-^ 
and our equation is 

(tr S+1){A — A(S+ l)" 1 - (tr S+ 1)~ 1 A\{S+ 1) = 0. 

The solution got from the first factor is clearly identical with that got from the 
third. Taking the remaining factor we have, on doubling it, 

A—2A(S+ l)- 1 + A — (tr / S+l)- 1 A = 

or {A-2A(S+l)- 1 \ + tr{A—2A(S+l)- 1 \ = 0, 

the solution of which is evidently 

A — 2A(S-\- l) _1 = any zero-axial skew matrix. 

Denoting by — N such a partially arbitrary matrix — that is to say, a matrix 
involving as was desired in(n— 1) different elements — we have 

2A(S+1)~ 1 = A + N, 
and .-. h(S+ l)A- 1 = (A + N)-\ 

and finally S= 2 (A + #)~ X A — 1 . 

7. We have thus reached the following theorem : 
The transfmmation of 

(x, y,z, \A\x, y,z, ) into (£,>', £, . . . .\A\%, n, 

can be effected by the substitution 

(x, y,z,... .) = (2(A + N)- 1 A- lJS, *, £ ), 
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where N is any zero-axial skew matrix arbitrarily chosen except for the condition 

that the determinant of the matrix A + N must not be zero. 

Or, following Cayley's example when dealing with the much simpler case of 

an orthogonal substitution, we may put the result in the form of a "rule," viz. 

To obtain the coefficients a rs of the linear substitution which transforms the 

quadric 

ax 2 + by % + cz 2 + . . . . + 2fyz + Igzx + 2hxy + 

into itself, form the determinant K adjugate to 
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or A 



and then put 
and 



when r =f= , 



A . a r5 = 2 (r th col of K\s tn col. of A) 

A . a„ = 2 (?■"' col. of K\r th col. of A) — A, 

A being the discriminant of the quadric. 

This is manifestly just as concise as Cayley's rule for the very simplest case. 
It is seen to pass into the latter when a = b = c= .... =1 and /= g = h 

— ■■ ■ • • • • — \J • 

8. It should be noticed that the equation 

trS.A.S=A 

resembles an ordinary algebraical equation which has only even powers of the 
unknown ; for, if S=(i be a solution, it is clear that S= — (i is also a solution. 
Indeed instead of resolving tr S .A . S — A into 

(trS+l)\A — A(S+ I)" 1 — (trS+l)- 1 A\{S+ 1), 

we might have resolved it into 

(tr S— 1) \A — A (1 — S)- 1 — (1 — tr S)-^ }{S— 1), 

the solution tr S= 1 being, to say the least, just as manifest as the solution 
tr S=—l. 

The roots obtained are thus : 

± 1, ±|2(^ + JV)-U— If. 
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9. Another point worthy of attention is the fact that the identity of 

(trS+l)\A — A(S+l)- 1 — (trS+iy'AKS+l) 
and trS.A. jS — A 

is quite independent of the form of A, and that therefore whatever A may be, 

the root of 

A — A(S+ l)- 1 — (fr/S'+l)- 1 ^ = 

is a root of 

trS.A.S—A = 0. 

When, however, we proceeded, towards the close of §6, to solve for the root of 
the former equation, we introduced the condition A = trA: consequently when 
A is unconditioned* it cannot be expected that the root thus reached will be a 

root of the equation 

trS.A.S—A = 0. 

10. Returning now to the general equation with which we started, viz. 

triS 2 .D.S 1 = D, 

let us give /Si a form analogous to the form 2 ( A + N)~*A — 1 obtained for S 
in §6, changing A, of course, into the more general D and, in order that we may 
have the full number of arbitrary constants, putting a perfectly arbitrary matrix 
M in place of the zero-axial skew matrix N. Our equation then becomes 

tr g 2 . D .{ 2 (D + M)~ l D — 1 } = D 
or trS % .\2D{D + ]\i)- l D — D\ =D 

or trS^lZDiD + irj-i-llD = Z>, 

and /. tr S i = \2D(D + M)' 1 - 1\~\ 

11. In place therefore of the simple theorem of §5 we have the following as 
a useful alternative : 

The transformation of the bipartite quadric 
(x, y,z, ... .IDfx', y', z' ) into (%,>!,%, W, yi\ £'» ) 

* The solution of A — -AOS+l) -1 — (tr S+1)~'A = 0, when A is unconditioned, does not seem to 
be easy. As we see readily from the complete equation tr S.A.S=:A, there are then implied n 2 
equations for the determination of n 2 unknowns. 
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ofr, as the saying is, into itself, can be effected by the substitutions 

{x,y, z, ) = (2A~ 1 D — 1#£, >?,£, ), 

(x\ y>, z',....) = tr (2D A- 1 - ljf , y,', % ), 

where A = D + an arbitrary matrix and has a non-vanishing determinant. 

1 2. Bearing in mind the equation tr S z . D . $ = D, we see that the foregoing 
result rests finally on the identity 

( 2DA~ 1 — 1)- ID ( 2A- l D — 1) = D , 

which is the same as 

Z)(2A" 1 Z>— 1) = (2DA- 1 — 1)Z>.* 

This latter however leads with equal naturalness to the identity 

(2DA- 1 — 1)D(2A~ 1 D—1)~ 1 = D, 

so that another form of the solution of the equation 

trS 2 .D.S 1 = D 



is available, viz. 



$ = (3A-ID-1)-M 
S 2 = tr(2DA~ 1 — 1)J 



13. The existence of an alternative form of solution is, of course, what 
might have been expected from the character of the function under discussion, 
which itself has two forms : 
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* The ultimate basis is, of course, the fact that 

WA-VD — D 



is equal to either 



I)®*- 1 !) — 1) or (32M- 1 — 1)2>. 
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Consequently, when we obtained the solution 

5 1 = 2A~ 1 D — 1, | 

5 2 = tr(2DA- 1 -l)- 1 ,} 

we might, by attending to the evidently permissible interchange, have foretold 

the solution 

S z =2trA- 1 trD—l, ) 

/Si = tr (2 trDtrA- 1 — l)~\ ) 

And this latter will be found to be identical with the second solution obtained 
above ; for — to take only the value for /Si — by the laws of matrix algebra 

tr (2 tr D tr A" 1 — l)~ 1 = \tr (2 trDtrA' 1 — l)}- 1 

— {tr (2 trDtrA- 1 )— l}" 1 
= (2A~ 1 D — l)- 1 . 

14. Further, however, it must be noted that there is no reason a priori for 
the unlikeness — which is notable in both of these solutions — between the value 
of /Si and the value of /S 2 ; and the presumption to which we are thus led that 
either form must be suitable for both, is fully borne out on investigation. Thus, 
taking the form for tr S. z in the first solution, we have 

(2DA- 1 — l)- 1 =\(2D — A)A- 1 \-\ 
= A{2D — A)-\ 
= \2D—(D — M)\(D — M)~\ 
= 2D(D — M)- 1 —1, 

which closely resembles the form for /Si in the same solution. Similarly, taking 
the form for /Si, we have 

2A~ 1 Z> — 1 = A" 1 (2D — M), 
= A- 1 (D — M), 
= \{D-M)- 1 A\-\ 
= {(D — M)- 1 {2D — D — M)\-\ 
= \2D(D — M)- 1 — \\~\ 

which is quite similar to the form for tr S. z in the same solution. 

15. And this is not all, for, since 

(D + M)- 1 (D + M) = 1, 
and .-. (D + M )~ X D + (D + M)~ 'M= 1 , 
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it follows that 

2(D + M)~ 1 D— 1 or S 1 = l — 2(D + M)- 1 3f, 

and similarly, 

2D (D — M)~ x — 1 or tr $ = 1 + 231 (D — 31)~ \ 



It also appears incidentally in the preceding paragraph that 

S 1 = (D + M)- 1 (D 
and irS z =(D + 31) (D- 



S 1 = {D + 31)- 1 (D-M) ,1 
i — M)- 1 .} 



16. We have therefore four identical expressions for $ in our first solution, 
and four identical expressions for $j, each of the former being quite similar in 
form to one of the latter, viz. 

S 1 =2(D + M)- 1 D — 1, and tr S,= 2D(D — M)~ l — 1, 
= 1 — 2{D + M)~ 1 M, =l + 2M(D—M)-\ 

= (£> + Jli")- 1 (Z) — 31), =(D + M){D — M)-\ 

= \2(D— M)~ l D—l\~\ =\2D(D + 31)~ 1 — 1 \~ K 

Looking at the first and last of the four forms we see that each of the eight 
expressions is such that to change the sign of M in it is the same as to take its 
reciprocal. Since therefore, in the second solution of §13, the values of $ and 
$j are respectively the reciprocals of S t and S 2 in the first solution, it follows 
that the second solution can be got from the first by merely changing the sign 
of M. But 31 is entirely arbitrary ; consequently the second solution is not 
really different from the first. 

17. The only advantage of these more complicated forms of S x and S% is 
that specialization from them is very easy. 

Thus, when D is axisymmetric, equal to A say, and Si = S 2 , we have for 
the conditioning of 31 

2 (A + M)~ l A — 1 = tr { 2A (A — M)~ l — 1 \, 
whence (A + 31 )~ 1 A = tr(A — 31)~ J tr A , 

and .-. (A + 31)' 1 = {tr A — tr 31)- \ 

and finally, 31= — tr 31, 

which is the definition of a zero-axial skew matrix. 
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18. In order that previous work on the above matters may be known, I 
append hereto a first approximation to a complete list of writings on Matrices. 
A supplementary list, to which I invite contributions, will be published when it 
is sufficiently bulky to warrant attention. 

LIST OF WRITINGS ON THE THEORY OP MATRICES. 

(1857-1893.) 



1857. 

Cayley, A. A memoir on the theory of matrices. Phil. Trans. R. S. London, 
CXLVIII, pp. 17-37; or Collected Math. Papers, II, pp. 475-496. 

Cayley, A. A memoir on the automorphic linear transformation of a bipartite 
quadric function. Phil. Trans. R. S. London, CXLV1I1, pp. 39-46 ; 
or Collected Math. Papers, II, pp. 497-505. 

1867. 

Laguerre, Sur le calcul des systemes lineaires. Journ. de l'Ecole Polj't., 

XXV, pp. 215-264. 

1870. 

Peirce, Benjamin. Linear associative algebra. American Journ. of Math., IV, 
pp. 97-215. 

Peirce, Chas. S. Description of a notation for the logic of relatives. Mem. 
Amer. Acad. Sci., IX. 

1872. 

Spottiswoode, W. Remarks on some recent generalizations of algebra. Proc. 
Lond. Math. Soc, IV, pp. 147-164. 

1877. 

Frobenius, Ueber lineare Substitutionen und bilineare Formen. Crelle's 

Journ., LXXXIV, pp. 1-63. 



1879. 



(a b) 



Cayley, A. On the matrix \ I and in connection therewith the function 

\c d\ 

ax + h Messenger of Math., IX, pp. 104-109. 
ex + a 

30 
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1882. 

Sylvester, J. J. Sur les puissances et les racines de substitutions lineaires. 
Oomptes Rendus . . . . , Paris, XCIV, pp. 55-59. 

Sylvester, J. J. Sur les racines des matrices unitaires. Comptes Rendus 

Paris, XCIV, pp. 396-399. 

Sylvester, J. J. On the properties of a split matrix. Johns Hopkins Univ. 
Circular, pp. 210, 211. 

1883. 

Buchheim, A. An identity in the theory of matrices. Messenger of Math., 
XIII, pp. 65, 66. 

Sylvester, J. J. On the equation to the secular inequalities in the planetary 
theory. Phil. Mag. 5th ser. XVI, pp. 267-269. 

Sylvester, J. J. On the involution and evolution of quaternions. Phil. Mag. 
5th ser. XVI, pp. 394-396. 

Forsyth, A. R. Proof of a theorem by Cayley in regard to matrices. Messen- 
ger of Math., XIII, pp. 139-142. 

1884. 

Sylvester, J. J. Lectures on the principles of universal algebra. Amer. Journ 
of Math., VI, pp. 270-286. 

Sylvester, J. J. On the solution of a class of equations in quaternions. Phil 
Mag. 5th ser. XVII, pp. 392-397. 

Sylvester, J. J. Sur les equations monoth6tiques. Comptes Rendus .... 
Paris, XCIX, pp. 13-15. 

Sylvester, J. J. Sur l'equation en matrices px-=.xq. Comptes Rendus .... 

Paris, pp. 67-71. 
Sylvester, J. J. Sur l'equation en matrices px = xq. Comptes Rendus .... 

Paris, XCIX, pp. 115, 116. 

Sylvester, J. J. Sur la solution du cas le plus general des Equations lineaires 
en quantites binaires, c'est-a-dire en quaternions ou en matrices du 
second ordre. Comptes Rendus . . . . , Paris, XCIX, pp. 117, 118. 

Sylvester, J. J. (Theorem regarding two matrices which have a latent root in 

common.) Educational Times, XXXVII, pp. 297, 385 ; or Math, from 

Educ. Times, XLII, p. 101. 
Sylvester, J. J. Sur la resolution generate de l'equation lineaire en matrices 

d'un ordre quelconque. Comptes Rendus . . . . , Paris, XCIX, pp. 

409-412. 
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Sylvester, J. J. Sur la resolution generale de liquation lineaire en matrices 
d'un ordre quelconque. Comptes Rendus, . . . . , Paris, XCIX, pp. 
432-436. 

Sylvester, J. J. Sur les deux methodes, celle de Hamilton et celle de l'auteur, 
pour resoudre liquation lineaire en quaternions. Comptes Rendus 
, Paris, XCIX, pp. 473-476. 

Syivester, J. J. Sur l'achevement de la nouvelle m6thode pour resoudre 
F equation lineaire la plus g6nerale en quaternions. Comptes Rendus 
, Paris, XCIX, pp. 502-505. 

Sylvester, J. J. Sur 1 'equation lineaire trinome en matrices d'un ordre quel- 
conque. Comptes Rendus . . . . , Paris, XCIX, pp. 527-529. 

Buchheim, A. Proof of Sylvester's " Third Law of Motion." Phil. Mag. 5th 
ser. XVIII, pp. 459-460. 

Sylvester, J. J. Sur la solution explicite de 1 'equation quadratique de Hamil- 
ton en quaternions ou en matrices du second ordre. Comptes Rendus 
, Paris, XCIX, pp. 555-558. 

Sylvester, J. J. Sur les conditions de l'existence de racines egales dans l'equa- 
tion du second degre de Hamilton, et sur une methode generale pour 
resoudre une equation unilaterale de n'importe quel degre en matrices 
d'un ordre quelconque. Comptes Rendus . . . . , Paris, XCIX, pp. 
621-631. 

Buchheim, A. On the theory of matrices. Proc. Lond. Math. Soc, XVI, pp. 
63-82. 

Cayley, A. On the quaternion equation qQ — Qq' = o. Messenger of Math., 
XIV, pp. 108-112. 

Sylvester, J. J. On Hamilton's quadratic equation and the general unilateral 
equation on matrices. Phil. Mag. 5th ser. XVIII, pp. 454-458. 

1885. 

Bochheim, A. A theorem in matrices. Messenger of Math., XIV, pp. 167, 168. 

Cayley, A. On the matrical equation qQ — Qq' = 0. Messenger of Math., XIV, 
pp. 176-178. 

Weyr, Ed. Sur la theorie des matrices. Comptes Rendus . . . . , Paris, C, pp. 

787-789. 

Muir, Thomas. New relations between bipartite functions and determinants, 
with a proof of Cayley's theorem in matrices. Proc. Lond. Math. Soc, 
XVI, pp. 276-286. 

"Weyr, Ed. Repartition des matrices en especes et formation de toutes especes. 
Comptes Rendus . . . . , Paris, C, pp. 966-969. 
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1886. 

Buchheim, A. An extension of a theorem of Professor Sylvester's relating to 
matrices. Phil. Mag. 5th ser. XXII, pp. 173, 174. 

1887. 

Buchheim, A. On a theorem of Professor Klein's relating to symmetric matrices. 
Messenger of Math., XVII, p. 79. 

Buchheim, A. Note on matrices in involution. Messenger of Math., XVIII, pp. 
102-104. 

1890. 

Taber, Henry. On the theory of matrices. American Journ. of Math., XII, 
pp. 337-396. 

Taber, Henry. On certain identities in the theory of matrices. American 
Journ. of Math., XIII, pp. 159-172. 

Taber, Henry. On the application to matrices of any order of the quaternion 
symbols S and V. Proc. Lond. Math. Soc, XXII, pp. 67-79. 

1891. 

Taber, Henry. On certain properties of symmetric, skew symmetric, and 
orthogonal matrices. Proc. Lond. Math. Soc, XXII, pp. 449-469. 

1892. 

Metzler, W. H. On the roots of matrices. Amer. Journ. of Math., XIV, pp. 
326-377. 

Metzler, W. H. On certain properties of symmetric, skew symmetric and 
orthogonal matrices. Amer. Journ. of Math., XV, pp. 274-282. 

1893. 

Taber, Henry. On the linear transformation between two quadrics. Proc. 
Lond. Math. Soc, XXIV, pp. 290-306. 

Metzler, W. H. Compound determinants. Amer. Journ. of Math. [XVI, pp. 
131-150], pp. 149-150. 

Taber, Henry. On orthogonal, substitutions that can be expressed as a func- 
tion of a single alternate (or skew symmetric) linear substitution. 
Amer. Journ. of Math., XVI, pp. 123-130. 



